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Abstract 

Explicit divergences and counterterms do not appear in the differential renormalization method, but they are 
concealed in the neglected surface terms in the formal partial integration procedure used. A systematic real space 
cutoff procedure for massless </)* theory is therefore studied in order to test the method and its compatibility with 
unitarity. Through 3-loop order, it is found that cutoff bare amplitudes are equal to the renormalized amplitudes 
previously obtained using the formal procedure plus singular terms which can be consistently cancelled by adding 
conventional counterterms to the Lagrangian. Renormalization group functions P{g) and j{g) obtained in the 
cutoff theory also agree with previous results. 

1. Introduction 

The purpose of this paper is to explore another aspect of the recently proposed differential renormalization 
procedure [1]. That method relies on the observation that (essentially all) primitively divergent Feynman graphs 
are well defined in real space for non-coincident points, but too singular at short distance to allow a Fourier 
transform. A regularization procedure must supply a prescription for the real space amplitudes which defines the 
short distance singularities such that integrals over them are well defined. This was done in [1] by a method 
which simultaneously regularizes and renormalizes amplitudes. The ideas involved are quite simple and best stated 
in terms of the 1-loop 4-point bubble graph of massless <f>* theory in 4 dimensions. This involves the singular 
function 1/x'^, where x'^ = {x^x^) , and is regulated as follows: 

1. Express such singular functions as derivatives of other functions which have well defined Fourier transforms. 

For example, 

1 1 InK'Px^ , , 

— = -7° ^- (11) 

X 4 x^ 

is an identity for x ^ Q, and the function In (M^x^) /x"^ has Fourier transform — 47r^ln yp^/M J /p^ where 
M = 2M/7 and 7 = 1.781 ... is Euler's constant. 

2. Use formal partial integration of the derivatives in (1.1) to compute integrals such as the Fourier transform. 
Thus the regulated Fourier transform of 1/x'^ is defined as — tt^ In ip^/M J . 

In Ref. [1] it was shown in a very explicit study of massless (f)^ theory through 3-loop order, that these 
ideas can be extended to renormalize all IP! vertex functions, including both primitively divergent graphs and 
those with divergent subgraphs. It was also shown that the resulting amplitudes satisfy the renormalization group 
equations in which M appears as the expected scale variable. Further applications of differential renormalization 
to gauge theories, supersymmetry, and amplitudes with massive particles have recently appeared [2]. 

Explicit divergences and the counterterms which cancel them never occur in differential renormalization. 
The usual ultraviolet divergences of field theory are hidden in the short distance surface terms which are dropped 
in step 2 above. It was an implicit article-of-faith in [1], justified only in 1-loop order, that these surface terms 
could be cancelled by counterterms for wave function, mass, and coupling renormalizations. Since this is crucial 
to the consistency of the procedure, we undertake to demonstrate it here up to 3-loop order in (p^ theory. For 



this purpose we wish to repeat the calculations of [1] using an explicit cutoff, implemented by modifying the 
Euclidean massless scalar propagator as follows, 

Real space calculations with this propagator are modelled as closely as possible on the differential methods of [1], 
leading to bare amplitudes r''{xi,e). In the limit of small e, we show that the bare amplitudes for each diagram 
can be expressed as the renormalized amplitudes of [1] V^ixjM) plus additional singular terms involving 1/e'^ 
or (ine^M^) . The latter are cancelled by adding local counterterms to the Lagrangian and including graphs 
generated by counterterm vertices. The scale M is required for dimensional reasons in the separation of regular 
and singular terms as e — > 0. It will be clear from our calculations that the singular terms are related to the 
surface terms neglected in [1], and the consistency of Step 2 above is thereby demonstrated. 

This investigation was primarily motivated by skeptics of the methods of [1] who were not convinced that 
overlap divergences were treated correctly and suspected an attendant violation of unitarity. We believe that the 
present investigation resolves such doubts at the concrete level of current calculations. Specifically, the fact that 
the singular cutoff dependence of bare amplitudes is cancelled by a Hermitean modification of the Lagrangian 
effectively proves that the results of [1] are unitary, provided that any non-unitarity in the cutoff chosen vanishes 
as the cutoff is removed. In the present case the Fourier transform of the cutoff propagator is 

and contains a logarithmic branch point at p = 0, with the cut corresponding to time-like Lorentzian momentum. 
This behavior is displayed in the limiting form for small e and fixed p, 



P P 



Ki{ep)-- + je' 



ln(^^^l-l 



0{e'). (1.4) 



Unitarity is satisfied as e ^ since effects of the logarithm vanish quadratically and the quadratic divergences in 
cutoff amplitudes are directly cancelled by additive mass counterterms, so there are no 1/e^ terms which multiply 
(1.4). 

One should note that -Ki{ep) falls exponentially as p ^ oo. Our real space computations are therefore 
equivalent to a momentum space approach with a damped propagator similar to that used in [3] to study the 
renormalization group in quantum field theory. 

Renormalization group equations are usually derived [4] by studying the relation between bare and renor- 
malized amplitudes. Since we have now systematically defined cutoff bare amplitudes within the differential 
renormalization method, we can repeat this derivation and obtain (3{g) and 7(5) as conventional scale derivatives 
of the renormalization constants. Our results agree with those found using the "experimental" approach to the 
renormalization group equations taken in [1]. 

2. The Cutoff Method in 1-Loop Order 

The Lagrangian of massive Euclidean signature cf)^ theory is 

C^\{d^f + \m^<i^^ + ^^\cpK (2.1) 

We will be concerned only with the massless case in this paper. The bare amplitudes for all Feynman diagrams 
will be computed using the Feynman rules of (2.1) with the cutoff propagator (1.2). These Feynman diagrams 
are shown in Figs 1-2, in which the same designation as in [1] is used. 



See figure 1 of Reference [1] 

Fig. 1: Diagrams which contribute to T^'^^x — y) in cf)^ theory. The divergence associated to the tadpole diagram 
b is immediately cancelled by an appropriate counterterm in our renormalization scheme, so 2- and 3- 
loop graphs which include tadpoles need not be considered. 

The classical contribution to the 4-point function, diagram e, and the bare amplitude for the 1-loop bubble 
graph f , are given by 

Teixi) = -167r^g(5i25i3(5i4 (2.2) 



and 



where 



and the notation 



V){x^,e) = 852512(534— 2 + 2-perms (2.3) 

[xfa + £2] 

Xij ~ Xi — Xj 

Oij — O yX{ Xj j 

will be used throughout. 

As in Section I of [1], we express l/{x'^ + e^)^ as nG{x'^,e'^), leading to the ordinary linear differential 
equation 



1 4 d f odG 



.2 



z 



(z + e^Y' z dz \ dz 
in the variable z = .t^. The general solution is 



(2.6) 



1 ln[(.T2 + e^)le'^\ a_ 
4 x^ x^ 



C(^2^,2)^__r:iv:__^^ + __ + ;,, (2.7) 

The additive constant b is irrelevant and can be dropped. The basic differential identity then becomes 



1 1 \n[{x^+e^)/e^] 

r° 15 ■ (^-^i 



(x2+e2)2 4 

We have chosen a = because the right hand side would otherwise have a S{x) singularity not present on the 
left. One also sees that the behavior of G [x^,e^) is sufficiently soft as a; ^ 0, so that the derivatives in (2.8) 
can be freely integrated by parts without generating a short distance surface term. 

We now introduce the dimensional constant M'^ and separate (2.8) into the two terms 

1 1 ln(x^+e^)M^ 1, 2,^2 1 /^^N 

- -D ^ ^> + -\ne^M^n-^. (2.9) 



(x2 + £2)2 4 a:;2 4 x 

Here D can be interpreted as operating to the right or left. Note that 

nX: = -4:Tr^S{x) . (2.10) 

As £ ^ we obtain 

^^ ^ -i □ i^ - ^ Ine-MHix) , (2.11) 

(x2 + £2)2 4 a;2 

in which we have introduced the notation Q to indicate that the derivative must now be interpreted as acting to 
the left in integrals, since ill defined singularities would otherwise be obtained. 



See figure 2 of Reference [1] 

Fig. 2: Diagrams which contribute to T^^' {xi, X2, x^, X4) in (f)"^ theory. 



To see the consistency of this interpretation, and its compatibility with Step 2 of the differential renormal- 
ization method let us compute the Fourier transform of (2.9) and compare with that of (2.11). The transform 
of the left side of (2.9) is easily evaluated using parametric differentiation, 



d^'xe'P-'' ^ == ^ / d^xe'P-'' 



1 9 /■ ,4 i„.^ 1 



(a;2+£2)2 2edeJ {x^ + e^) 

__27r2_a / 
s de \ 



-^1 (^^^^^^V (2.12) 



= 7r2p2 J [Ko(ep) + \neM] ^- 

^pZ pZ 

We have taken the simple exact result in the third line, and rewritten it with a mass scale M introduced, so that 
the final form corresponds to the Fourier transform of the two terms on the right side of (2.9) with D interpreted 
as the factor —p^. The limiting form of this result as e ^ 0, namely, 

-TT^lnp^/M^ -TT^lne^M^ (2.13) 

coincides with the Fourier transform of the right side of (2.11). 

The physical interpretation of these manipulations can be seen by substitution of (2.11) in (2.3) which then 

reads, 

• In x'^ M"^ 
T){xi,£) — > -2.9^(512(534 □ —^ h 2-perms 

"^13 

-24^V<5i2^i3<5i4ln£2M2 (2.14) 

= r^(x,,M) - 247rV'5i2(5i3<5i4lne^M2 . 

This is the sum of the renormalized amplitude of [1] plus a local term with singular coefficient Ine^A/^ which is 
removed by adding the counterterm 

-TrVlne^M^^* (2.15) 

to the Lagrangian (2.1). One can also use (2.13) to obtain the same interpretation in momentum space, and the 
renormalized amplitude, obtained from the first term of (2.13), agrees with the formal partial integration rule of 
Step 2. 

The other 1-loop diagram is the tadpole contribution h to the 2-point function. Using the damped propagator 
(1.2) one finds the bare contribution to the IP! 2-point function 

Tl{x-y,e)^%S{x^y). (2.16) 

This can be cancelled immediately by adding the mass counterterm 

-4^2 (2.17) 

to (2.1). Insertions of the tadpole in other diagrams are also cancelled by the counterterm graph from (2.17), 
and higher order tadpoles can be cancelled similarly. For this reason we do not consider tadpole diagrams further 
in this paper. 



The main purpose of this section was to show in 1-loop order using the cutoff procedure that the bare 
amplitude can be expressed for small e as the renormalized amplitude of [1] plus a singular term in e which can 
be absorbed by coupling renormalization, and that this term corresponds closely to the singular surface term 
neglected in the formal partial integration prescription of differential renormalization. 

3. 2-Loop Amplitudes 

We now apply the cutoff method of calculation to 2-loop diagrams, where we shall encounter some new 
features. 

Graph c is the first non-trivial contribution to the 2-point function. Its bare amplitude is 

'•'<-•'"' -IS (i^^ f^-') 

To handle this we simply differentiate the identity (2.8), obtaining 

1 1 ln(a;2+e2)/£2 32 

- DD^— ^^^ + -^-5-.. (3.2) 



(x2+e2)3 32^^ a;2 (a;2 + £2)4 

The second term in (3.2) is a function whose limit as e ^ vanishes for x ^ 0, but is singular for e = 0. One 
suspects that its limiting form is a distribution, and this may be confirmed by studying the integral 

d^xf{x)ie^ 3 fd^yfiey) 



(x2+£2)4 ^2 J (y2 + l)4 ' 

where f(x) is a smooth function which is damped at large distances. After the change of variables x — ey, one 
can evaluate the integral for small e by expanding f{ey) in a Taylor series through second order, and noticing 
that the contribution of the remainder is a convergent integral at large distance which vanishes as £ ^ 0. After 
explicit evaluation of two elementary integrals one finds the result 

which is equivalent to the statement 

l^, (^2f72)4 - ^S(-) + Y°^i-) ■ (3-5) 

Such "representations of distributions" appear frequently in our work, and they are collected systematically in 
the Appendix. 

The final step in the treatment of (3.2) is to introduce the scale M in the first term, obtaining 

nn'-^^^^^^P^^u6'^^+A.Hne^M^USix) (3.6) 

a;^ X'' 

by steps similar to those leading from (2.8) to (2.11). We now combine (3.4) and (3.6) and insert the result in 
(3.1) obtaining 

r^(x, e) - J^ DD i^ + f^ (insHl^ l) nS(x) ^^5{x) . (3.7) 

The first term is the renormalized amplitude r^(x,M) of [1] while the last two terms can be cancelled by wave 
function and mass counterterms. 

We now turn to the 4-point function, for which graphs g and h contribute in 2-loop order. The bare amplitude 
for graph g is 

b 4g^ 

r (xi,£) = — V'5i2'^344(a;i3) + 2-perms , (3.8) 



where /e(x) is the convolution integral evaluated exactly in (A. 2). The limiting form for small e is given by (A. 4), 
and we rewrite it as 



Ie{x) 



7r2 . lii'^x^M^ 
^D ^ 



Ine^M^ 



lna;2Af2 



+ Att^ In s^M^S{x) 



tt" In^ e^M'5{x) 



(3.9) 



When inserted in (3.8) the first term gives the renormalized amplitude F^, and the second term is proportional 
to the limiting form of r^(a;i,e) in (2.12) times the singular coefficient Ine^M^. This non-local divergence will 
be cancelled, as we show systematically in Section 5, by the bubble graph generated by the coupling counterterm 
(2.15). The last term in (3.9) can be cancelled by an order g^lv? e'^M'^4>'^ counterterm in the Lagrangian. The 
amplitude obtained after insertion of (3.9) in (3.8) is presented in Table 2. 

Graph h requires a longer discussion because it is the first diagram whose basic form is triangular. The bare 

amplitude is 

8.g3 1 1 1 



T,{x,,e) ^ -r^J,,-^—-, -^—-^ ^^^^ ^ ^^^^ + 5-pcrms 



(3.10) 



We let .xi4 = X, X34 = y. Using (2.9) or (A. 10), and (2.10) together with the antisymmetric derivative 
identity 

ADB = df,(Ad^B\+BnA, (3.11) 

one sees that 



{x — yY + £^ (a;^ + e^) (y^ + e^)^ 

1 
- tD 



1 d 

49^ 



d ln(j/2 + e2)M2 



[x — yY + £^ x^ + £^ dy^ 



ln(j/^ 



e')M' 



4 (X - yY + £2 (a;2 -I- £2 



7r2j(y)lnMV 



y 



1 lnM^(x^ + £^) 



+ 7r2lnMV(5(a;) 



(3.12) 



Note that (2.9) was used both to replace l/iy^ + e^Y on the left side of (3.12) and to replace l/(a;^ + e^Y in 

the last term. The first term is regular as e ^ provided we understand that -^ must be integrated by parts. 
A similar remark applies to In Af^(a:;^ + e^) term in the last line. 

The second term in (3.12) is an example of something we call a triangular structure. To study its limit as 
£ ^ 0, we use the simple identities 

1 =.2 



1 



{x — yY + £^ {x — yY {x — yY [{x — yY + £^ 



(3.13) 



= -47r^(5(x- y) -D- 



(x — yY + £^ [x — yY [{x — yY + £^] 

Using (3.14) the second term in (3.12) can be written as 

£2 In (y2 + £2) M2 



IT S{x — y) 



In(x2+£2)M2 1 



D 



(x2+£2) 4 (x-y)2((a;-y)2 +£2) (.x2 + £2)y2 



(3.14) 



(3.15) 



The limiting form of the first term is obtained from (A. 11a) which is essentially a differential identity. 

We claim that the limiting form of the second term is that of an eight dimensional delta function 
— jCS{x)S{y) corresponding, after insertion in (3.10), to another two-loop order coupling constant countert- 
erm. This claim can be verified by studying the integral of the term in question with a test function f{x,y). 



After scaling variables, x -^ ex, y -^ ey, one sees that the limiting contribution involves only /(O, 0), and that 
the constant C is given by 

= -Att^ (In e^M^-B) 

where the first integral becomes trivial after partial integration of Dj,. The evaluation of the second integral is 
discussed briefly in the Appendix, see (A. 13), and the result of the argument beginning with (3.13) is given in 
(A. 12b). The identities (3.13) and (3.14) are useful because the product terms fall off fast enough in the infrared 
so that their contribution can be obtained by scaling arguments. 

It is well worth noting that when scaling arguments are used in the study of integrals involving test functions, 
the question of the limit as e ^ is effectively transferred to the question of the behavior of the large x, y behavior 
of integrals over the scaled variables. It is a correct rule of thumb, which can be verified by more careful limiting 
arguments, that the limiting contribution of a term in the bare amplitude is a J-function or product of (5-functions, 
if the integral determining the naive coefficient of the (5-functions is infrared convergent as is the case for B and 
C in (3.16). 

The results (A. 11a) and (A. 12b) are now combined with the simple limits of the first and third terms of 
(3.12) to obtain the limiting form of F^ given in Table 2. Again one finds the renormalized amplitude of [1] plus 
singular terms to be cancelled by counterterms. 

4. 3-Loop Diagrams 

We now continue the program of the last two sections and study the limit as e ^ of cutoff amplitudes for 
the 3-loop graphs shown in Figs. 1 and 2. It is worth emphasizing that bare amplitudes are independent of the 
mass scale M, 

M^r\x,,e)^0, (4.1) 

because M is introduced only to separate r''{xi,e) into regular and singular terms. Thus one can use the property 
(4.1) as a check on the intermediate steps of the calculation of a complicated amplitude. The same mass scale 
is used in all diagrams in order to agree with the renormalization scheme of [1]. 

Graph d is the only contribution to the 2-point function at 3-loop order. Its bare amplitude is 

r;j(x,e) = ^ „ "^ . f d^u ^ ^ ^ „^^ i 2 . (4.2) 

We now use the result (A. 2) for the convolution, except that we compute D acting on (A. 2). We obtain. 



2.g3 1 j {x^+2e^)\n{x^ + 2e^-\x\y/x^TIP) /2e'^ 

TT^x^+e^[ \x\3^x^+Ae^^ ' x"^ {x"^ + 4e^) 



T',{x,s)^-^-^-^{'- \ .. , . ' =3 ^^ + Z^TZ^-rST ^ • (4-3) 



Inspired by the form of the renormalized amplitude, we compare (4.3) and 

an — ^ ^ — ^ — 



327r2 
2g3 r(a;2-2e2)ln(.T2 + e2)/e 



7r2(a;2 + e2) !_ (x2+e2)3 ' (2.2+£2)3j • ^^-^^ 

The difference between (4.3) and (4.4) again involves a representation of the distributions S{x) and nS{x), as 
one can verify by integration with a smooth f{x) as in (3.3-3.4). One can then write 



rda^,e H-Trf^nn^ ^^ 5 ^ ^^ I^D.Six)- f^D2n6ix)+0{e) , 4.5 

where Di and D2 are the purely numerical values of the following integrals, 



, . 1 I (.t2 + 2) In (x^ + 2- \x\Vx^ + 4) /2 1 

Di = ax- ' ' 



1 i |a;|3ViM:4 a;2(x2+4) 

(,x2-2)ln(a;2 + l) x^ 



(,x2 + l)3 (,^2 + 1)3 J ' 

T \ ( X 

Do= I d^x 



,4 x2 j (x2 + 2)ln(x2 + 2-|a;|Vi4T4)/2 



-1 [ |a;|3V^M^^ a;2(2.2+4) 

(x2-2)ln(x2 + l) x2 ] 

(x2 + l)3 (,^2 + 1)3 J ■ y ■ ) 

The logarithmic terms in (4.5) may now be separated into regular and singular terms in e after introduction of 
the mass scale M. One then finds the limiting form given in Table 1. In this expression Di and D2 are the 
coefficients of finite counterterms at 3-loop order. These terms become relevant to the cancellation of divergences 
only at the 4-loop level, so the integrals (4.6) need not be evaluated. 



Table 1. Bare cutoff amplitudes for graphs contributing to the 1 PI 2-point function. The subscripts denote 
the graphs shown in Fig. 1, and the tadpole graph b is omitted for reasons discussed in Sec. 2. The first 
term in each entry is the renormalized amplitude obtained in [1], and this is followed by cutoff dependent 
terms in the limit of small e. The numerical constants Di and D2 are the values of the integrals in (4.6). 



Ta{x) = -n5\x) 



r!!(.,s) ^ J^ in i^ - f! (1 - lne^M2)n^^(.) - ^^^^(x) 



u, , g3 . . ln2.x2M2 + 31nx2Af2 

+3fflne2M2r,(.T,e) - ^ (in^ sHl^ + \n e'' M^ + '^\ n5\x) 



We now begin our treatment of the 8 graphs which contribute to the 4-point function. We shall be rather 
brief in our discussion of the easier graphs and concentrate on the more difficult ones, namely, j, 1, n, and o. 

Graph i is particularly easy since the bare amplitude 

i^ 4 /* 1 

rj(xi,e) = ^-r<5i2(534 / d'^ud'^v ^ 5 5- + 2-perms , (4.7) 

^^ J [{xi - uf + e-^f [{u - vf + e-^f [{v - x^f + e-^f 



is a double convolution of factors for which the identity (2.8) may be used. As discussed in the Appendix, the 
limiting form of the convolution is correctly given by the convolution of the limiting form (2.11) of each factor. 
This leads to the result given in Table 2. 

Graph j is somewhat more involved. We start by writing its bare amplitude, 



r''(xi,e) = H — T^i'i—n ^^5 ^ / d'^u 7T 7T + 5-pcrms . (4 

We substitute the result of the integral (A. 2) and then use (3.11), which splits the amplitude into two parts, 



^4 1 /I 



--^S,2^—-^03 -r— 72 ^ -2 + 5-pcrins . (4.9) 

We take the limit e ^ in the first one, and introduce M . This yields 



'f^'^''^'Ma:?,+e2:,2 ^ .2 ^3m 



In^ (xl^ + 2e2 _ |a;34iy^l7T4i^) /2< 



13 ^ ^ -^14 I ^ -^34 

Tt \'''13-^14 •''34 / ^ \-^13 •'•14 •'•34 

-fi^-^^^ (^^V^^^a.^^) We^M^ . (4.10) 

The first term of the R.H.S. is a piece of the renormalized amplitude of Ref. [1]. The second is a piece of 
g\ne^M^Tl{x^,e). 

The third term requires a little more work. We undo relation (3.11) recognizing a delta term plus the 
triangular structure (A. 12a), 

^7<5i293u ( .^ . . ^ . ^3.4- ] In^eHl^ 



TT^ ^ \ a;?, + e^ 2;2 j_ ^2 t^ ^2 



13 I = X]^4 -r t .^,34 



1 In^p^M^ 

4g%2'?34 . , ,. In' e'M' - 4/ <5i2(5i3 + 4g\Hi25i3Su In e^M^ . (4.11) 



Let us now turn our attention to the second term of (4.9), which is a triangular structure independent of 
scale M. It will have a representation of the following form in the limit as e ^ 0, 



5-(5i2n3 



TT^ \ xu + e^ J Xta + e 



g\ / 1 \ 1 ln[xl^ + 2e^-\x3i\^4^T4^)/2e^ 



13 -r c / j.]^4 -r c j.34 



.4;C._ J X _. /(^34) 

^34 (2^34 + ^ ) 



g''(5i2 <^ S 



Ji3- 



+ ^25j,5i3<5i4| , (4.12) 



where Bj is a numerical constant. We shall determine the function /(a;34) by the trick of comparing the left side of 
(4.12) with a similar expression in which the logarithm with complicated argument is replaced by In (x^ +e^)/e^. 

It can be checked that the difference between the two expressions amounts, in the small e limit, to a delta 
function whose constant coefficient Cj is given by 



In' (.t2 + 2e2 - \x\Vx^ + 4e^) /2e'^ - ln'(a;2 + e^)/e^ 



tt'c, = / d^x "' '-^ ■ -' — '^";,,:r,;;r — '"'^ ' ' "' ■ (4-13) 

9 



Thus, we have 



''34 



-%2^i3 "'rr2' ;;;;; + ^^'(s, + c,),5i2(5i3'5i4 . (4.14) 



^34(^34 + ^ ) 



We introduce M inside the logarithm, 



^4 



<5i2<5i 



^34(2^34 + ^ ) 



ln2(a;24 + £2)^2 _ 2 1n(x§4 + £^)M2 Ine^M ^ + In^ e2M2 
^34(2^34 + s ) 



4„ e ill V-^34~rc )iy± - z. 111^.^34 -p c ;iv>^ 111c ivi ^111 c ji-j 

4g 012013 2~r^ — n2^ • (4l3j 



Notice that the second term of the R.H.S. is another piece of glne^M^r^(a;i,e). The third term cancels with 
the second term of (4.11). Finally, we use (A. lib) to regularize the first term. These results are collected in the 
amplitude of Table 2, where bj = Cj + Bj. 

We now turn to graph k, whose bare amplitude reads 

Tf.{xi, s) = 

'^9^ . f A 1 1 / ^ 

— tOi2 I a u 5 • —^ — + 5— perms . (4.16) 

^^ J [(a;i - 7i)2 + e2]2 [(^.g _ y)2 + ^2] [(^^ _ y)2 + £2] [a;24 + e2] 

This is a convolution of bare amplitudes for the bubble and ice-cream cone subgraphs. However, the renormalized 
amplitude [1] for this graph is also a convolution of the corresponding renormalized amplitudes interpreted using 
formal partial integration. Since our earlier results show that the limiting forms of the bare amplitudes for graphs 
f and h are equal to the renormalized amplitudes with D or 9 derivatives plus singular terms, these expressions 
from Table 2 can simply be inserted in (4.16). One finds the renormalized convolution amplitude plus terms 
containing 5'^{u — xi) which render the ct^u integrals trivial, and this leads immediately to the result in Table 2. 

We now study graph (. starting from the bare amplitude 

Q 

^\{xi,e) ^ ^/(5i2'534i(a;i3) + 2-pcrms (4.17) 



i/here 



^(^) = 2 I 2 -^(^) 

x'^ + e 



I{x) = / cf^ud^- "^ ^ 



1 



n3 



[u — vY + e^ 



(4.18) 



[x — uY + S^ v"^ + £2 

The new problem that arises here is that the integral I(x) is infrared divergent, as may be seen by fixing u — v 
and considering the integral over u + v. The physical reason for this is that the subgraph c contains a mass 
shift which causes the subsequent integration over the massless propagators to diverge. One thus expects an 
infrared convergent result only when the mass shift counterterm insertion from (3.5) or (3.7) is subtracted. Our 
treatment will make this clear. 

We insert the identity (3.2) in (4.18), and study separately the two integrals 

1 /■ A A 1 1 In |"(u - 'y)2 + £2] /£2 

^^(-) - -32y-^^--^^- (.-.)2+.2.2+,2 °° ^' (/.)2 ^^ (4.19) 
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h {x) = -i ! d\ dt^v \. ^ -J-^ 4 • (4.20) 

The first convolution integral is infrared finite, essentially because the Fourier transform of the DD log term con- 
tains a factor of p^ which amply compensates for the l/p^ factors in the transforms of the other two propagators. 
Therefore we treat l2{x) first. 

It is not difficult to verify the following differential identity 

3e^ 1 e^ T^^ ,, X /.O.N 
1=0° ^ + 7r^Kz) (4.21) 

(z2+£2)4 g ^2(22 4.^2)2 2e^ ' ' 

in which the 5{z) singularity cancels between the two terms. We insert (4.21) with argument z ^ u — v \n 
(4.20). The contribution of the first term, called l2i{x), is infrared finite, again because there is a factor p^ jp 
momentum space from the D in (4.21). The infrared divergence is now isolated in the contribution 122(2;) of 
the last 5{u — v) term in (4.21), and it is clear that this integral will be cancelled completely when the mass 
counterterm for subgraph c is inserted. See (3.5). 

After partial integration of Di, and use of (A. 6), we find that l2i{x) can be written as 



hi{x) 



d^u d'^v e^ 1 

{X - UY + £2 (^;2 + £2)3 (^ _ ^)2 [(y _ „)2 + ^2^ 

d*ud*v 1 1 

{X - ew)2 + £2 ^y2 + 1)3 (y _ ^)2 [(y _ ^)2 _^ if 



(4.22) 



where we have scaled u ^ eu and v ^ ev \r\ the last line. It is legitimate to take the £ ^ limit inside the 
integral because the residual integral is finite. This gives the simple result 

/2i(.x)--^ (4.23) 



'ith 

C 



d'^u d*v 



{v^ + i)^{u-vy[{u-vy + if 



d'^z 



[v^ + iyj Z2(z2 + 1)2 



(4.24) 



= — TT . 

2 
We are really interested in the contribution of l2i{x) to L{x) in (4.18) which is given by the product 

^21(2;) = -^I^hiix) (4.25) 

and it is not correct to say that the limiting form of this product is obtained simply by inserting the limiting form 
of (4.23) of l2i{x). Instead we note the following general structure of l2i{x), namely 

/2i(x) = -V1f(£Vx2) (4.26) 

which follows simply from (4.22). Because of the result (4.23), we know that 

lim F {e^/x^) = 1 . (4.27) 

We can therefore write 
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which is an exact representation. Using (A. 9. a) one can see that the limiting form of the first term is that of the 
bare amplitude for the bubble graph f plus a S{x) term. Using again a test function and scaling argument, one 
can show that the limiting contribution of the second term in (4.20) is also of the local form C'd{x) where C" is 
a numerical constant defined by the infrared convergent integral 

A more explicit form can be found using (4.22), but is not necessary. 

The integral Ii(x) remains to be studied. We use the identity (3.13) for each propagator factor obtaining a 
representation with four terms 

1 s^ e^ 



h{x)= Jd^ 



'.d\ 



{x — u)^v'^ (x — uY [{x — u)^ + e^] v^ (x — u^-v^ (v^ + e^) 

{u — v) 



4 1 In [(U - W)2 + £21 /£2 

DD ^ , ,2 • (4-30) 



(x — m)2 [[x — m)2 + £2] y2 (-^2 _|_ ^2-j 

After partial integration ofOuOv the first term trivially becomes 

Ini.) = l6y^^^^^^^fi^. (4.31) 

Each of the last 3 terms in (4.30) has the structure -^F (^e^/x^) because of dimensional considerations. If 
we can show that F (^e^/x^) vanishes as x^ -^ oo, then the contribution of these terms to L{x) of (4.18) can 
be shown to be purely local by a scaling argument similar to that used for the second term of (4.28). We next 
discuss how to establish that F (^e^ /x^) vanishes for each of the last 3 terms. 

After partial integration of Dt,, the second term can be written as 



2 2/ d'^u ln[(u2+£2)/£2 



'-(")-^"^ / (x-.)2[(x-.)2+£2^ ° 



-167r2£2 



(4.32) 



(a; — m)2 [(a; — u)^ + e'^] (u'^ + £2) 



where (2.8) has been used. A crude estimate of the asymptotic behavior gives 

£2 f d'^U 

x^ J (u^ + e^y 

This is not quite correct since the residual M-integral is logarithmically infrared divergent, but it means that the 
correct falloff is 

p2 

Ii2ix)r^ ^Inx^/s^ . (4.33) 



The function F {^e^ /e^^ thus falls nearly a full power of e^ jx^ faster than necessary, so we are content with the 
heuristic argument above. The third integral in (4.30), namely /i3(a;), can be shown to be equal to Inix) after 
partial integration of Dm and change of variables, so I\?,(x) also satisfies (4.33). 

The fourth integral /i4(x) is more complicated, but it is not difficult to show that it falls rapidly as x ^ oo. 
If we crudely extract the factor l/x** and use (2.8) we find 

e^ /■ .^ ,A _ 1 1 



£"* /■ 



' [(m - W)2 + £2]2 i;2 („2 _^ ^2) 

^ / d^ud^vXIu 2 • ^v \^rr-o ^ ) (4.34) 

x^ ] [(w-i;)2+£2]2 V^'" («' + £')' 

4 /dVV i 2- fd^vv(^-4 ^ 

x^ J (u'2 + £2)2 J \v^(v2 + e^) 
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The translation of variables u' = u — v \s permitted because the eight dimensional integral is convergent. Each of 
the two four-dimensional convergent integrals in the last line vanishes by symmetry, and this means that Ii4{x) 
actually falls faster than l/x^. This is more than enough to conclude that its contribution to L{x) is purely local 
as e ^ 0. 

The total contribution of Ii{x) to L{x) is therefore 

L,{x) = 16^^ \, l^/f + C"S{x) (4.35) 

X'' (x^ + e^j 

where the first term comes from (4.31). We introduce the scale M in the logarithm and apply the differential 
identity (A. 11a). Similarly (A. 9a) and (A. 10a) are applied to the first term of (4.28). The limiting form of the 
bare amplitude r''g{xi,e) given in Table 2 is the contribution of these two terms plus the infrared divergent mass 
counterterm integral l22{x) discussed below (4.21) and a local triple 6-term. 

The bare amplitude of graph m is 

^ 4^4 1 1 11 

r„(x.,e) - —^-^-—^-^-^-—^^-—^-^—-^ + 5-perms . (4.36) 

We now use (2.9) for the bubble subgraph amplitudes which leads to 

., , / in(a;22+e2)M2 ln{xl4+e^)M^ 1 1 

r™(a^O = fuD , □ ^ o 9 , o ,2 , .9 + 5-perms 



"12 -^34 -^13 ^ ^ -^24 

_2 71 j-2-nb ^.^ „\ „2 i„2 „2 i\#2-pfc / 



+gliie^M^Tl{x„e) - g^ln e^M'^T'}{xi,e) . (4.37) 

Using (3.11) we obtain the antisymmetric derivative terms in Table 2, plus the expression 

47r yXj^g + e / x-^2 x^^^ £ a;^^ 

The limiting behavior of F can be obtained by a procedure involving the use of the identity (3.14) in two 
factors of (4.38), namely those involving D of the cutoff propagators with arguments 2:13 and 0:14. A detailed 
discussion would be lengthy, and since there are no essentially new techniques involved, we give only a brief 
description. It is convenient to split the product of logarithms in (4.38) into terms proportional to In e^M"^, 
Ine^M^, and Af-independent terms, noting that factors such as (in [(a:;2 + £2)/£2j^ ^^2 ^^^ non-singular as 
x ^ 0. Each of the three terms above gives rise to four terms coming from the product of the two identities 
based on (3.14), and each term can be studied separately in a straightforward way. A minor difficulty occurs 
in the ln^e27\,//2 term, because the product 1/ {xf2X'^4) becomes ultraviolet singular when the arguments are 
identified in the S(xis)S(x24) term of the identities (3.14). to handle this one uses essentially (3.13), namely 

1 1 £2 



•^12 -^12 T ^ '^12 V'''12 r ^ i 

The first term cuts off the ultraviolet singularity of the product 1/ {x\2x\^ . The contribution of the second term 
to the In e^IvP term of (4.38) is then studied, without use of (3.14), and can be shown to be of the local form 
(5i2i5i3(5i4 as e ^ 0. 

The result of the analysis above is the following formula for the e ^ limit of F: 

F ^ 4/<5i3<524i^4^^^2^^^ - 4/^2<5i3'5i2'5i4 (in^e^M^ - 2B\ne^M^ + S„) . (4.39) 

All triple (5-terms were obtained by studying the limiting behavior of integrals with a test function h (2:12, 2^13, 2:^14) 
of the three independent variables in (4.38). The purely numerical coefficient B^ comes from a combination of 
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integrals from various terms in the analysis above. One would expect that the singularity associated with the 
bubble subgraphs of m should involve the amplitude of the ice-cream cone subgraph h. It is therefore important 
to verify, as we have done in detail, that the coefficient B in (4.39) is given by the same integral, see (3.16) or 
(A. 13), that appeared in the original analysis of graph h. Only then we will have a consistent cancellation of 
divergences as e ^ by counterterms. A further check on (4.39) can be obtained by verifying that the limiting 
form of r^ given in Table 2 satisfies (4.1). The net coefficient of Blne^M^ is easily seen to cancel. 

Graph n is moderately complicated. Its bare amplitude is 



85 



4 



T':,{x,,s)^^S- 
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TT^ X?. + £2 2;2 . 



d*u ^ + 11 -perms . (4.40) 

[{Xi - ■u)2 + £2] [(2;^ _ u)2 _^_ £2] [(^.g _ y)2 ^ £2]^ 

We focus attention on the integral, use (2.8), and follow the spirit of the corresponding steps in [1] to obtain the 
two terms 

r^(x„ e) = -^<5l2^2^ ^ ^2 ^2^ + £2°3 J ^'" («2 + £2) ((^^^ _ y)2 + ^2) (^^3 _ y)2 ' 

In 71 ..^2 , ^2 +Tr'^i2 / 2 , ^2x/^2 , ^2x /^2 , ^2x2 + H-PCTms (4-41) 



Xxi4-?i)2+e2y 7^2 (a;23+£2)(2.2^+£2) (2,2^ ^ ^2) 2 

as an exact result. We now use (3.11) to split the first term above into two parts. The term with the anti- 
symmetric derivative is already finite as e ^ 0, so we can use the techniques of [1] to obtain the corresponding 
term in the renormalized amplitude. 

The second term resulting from use of (3.11) in (4.41) contains the difference of the two triangular structures 
(A. 12c) for n = 1 and (A.12d). The w-integral cancels and one finds that the total contribution of the first term 
in (4.41) has the limiting form 

-^Si2 k^ ( ^^ 53p In ^X(xi3, a;i4) ) + 11-pcrms + g^AB + A)T'f{x„e) (4.42) 

77 \2;i4a;34 Xi^ J 

where 

"'^^^y^-I'^'^ u^ix-uny-uy ^'-''^ 

is the same function introduced in [1]. K{x,y) has an ultraviolet finite Fourier transform. 

The remaining task is to study the second term in (4.41). This is straightforward if we use (A. 10c) and 
introduce the scale M to obtain the exact result 

g^ , 1 In2(x2 +e2)^^2^21n(x2 +£2)^2 
-hSi2T^-^Kr^-^°3 -^ 2 -^ + 11-pcrms 

+2g\ne^M^rl{x„e) - 2g'^ {\n^ e^M^ + 2lne'^M'^)T''f{xi,e) . (4.44) 

We use (3.11) again, obtaining an ultraviolet finite antisymmetric derivative term and a triangular structure, 
which is a combination of (A. 12. a - A.12.c). We then use (A. 11a) and (A. lib) and assemble our results to 
complete the amplitude given in Table 2. 

We shall use a different technique to analyze graph o in order to avoid an intractable interplay of internal 
integrals and e ^ limits. Namely, we will obtain the renormalized amplitudes plus non-local divergent terms 
directly, but we will use (4.1) to obtain the coefficients of scale dependent local triple 5 terms. This is a 
mathematically correct shortcut because the bare amplitude is independent of M. These considerations determine 
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the limiting form of amplitude except for a purely numerical multiple of (5i2<^i3<5i4 which is simply a change in 
the renormalization scheme of [1]. 

The bare amplitude of graph o is 



/o(^) 



d'^ud^v 



((u _ „)2 _^ g2)2 y2 _^_^2y2 _^_ £2 (3, _ y)2 + ^2 (^J - 77)2 ^ ^2 

We regulate the bubble subgraph factor using (2.9) and integrate Dm by parts obtaining the three integrals 



(4.45) 



/o(^) 



1 



1 



1 



47 \ u^ + e'^ v'^ + e"^ {x — uY + e"^ {x — v^ + e^ 



1 



1 



1 



1 



-2a, 



1 



y2 _|_ ^2 ^2 _|_ £2 (^2; — m)2 -^ e^ [x — w)2 + S^ 

1 „ 1 1 1 In ((u - 'y)2 + £2^ ]\^2 



A* 2 



:2^M 



y2 _|_ g2 ^2 _|_ g2 ^2; — m)2 + £2 (^2; — u)^ 



(w — vf- 



(4.46) 



plus a contribution to V\{xi,e) proportional to the bare amplitude of graph g, namely gXne^M'^V^Axi^e). 

The first two terms of (4.46) have the same form, containing the triangular structure (A. 12b). We are 
entitled to use formula (A. 12b) inside the integral provided that a; 7^ 0. To account for a possible singularity in 
this limit, we include a 8{x) term, with unknown coefficient F^^ieM). Such coefficient will be fixed at the end of 
the computation, requiring the amplitude to satisfy equation (4.1). So we have 



47r-^ 




1 



1 



1 



ln((u-i;)2+e2)M2 



112 + £2 (2; _ y)2 -^ ^"^ (x — u)2 + £2 

ln(^;2 + £2)m2 



(u — w)2 



+ 47r^(ln£^M2-B) 



1 



Fl{eM)b{x) 



(4.47) 



a;2 + £2 (x — u)2 + £2 il2(t|2 _)_ ^2^ J2;2 -|- £2]2 

We regulate the divergent term in the integral using (A. 11a), integrate by parts and regulate again using (A. 11a) 
and (A. lib). One thus obtains the following limiting form of the first two integrals in (4.46) 



In^ a;2 Af 2 ^ g In a;2M2 + 12 In x^M"^ 
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+7r* (- In^ e^M'^ - 2 In £^^"2 + 2S + 2) • 



(x2 +£2)2 



Fl{eM)6{x) 



(4.48) 



We now study the last integral in (4.45). If a; 7^ 0, the integral actually converges if the limit £ ^ is taken 

in the integrand, and the result 

lna;2M2 , 1 ' 



^■k" 



V 



(4.49) 



was obtained for this limit in [1] by mathematically correct steps not requiring formal partial integration. The 
role of the cutoff £ is therefore to determine the 8{x) term in the result of the integral. We are therefore entitled 
to assume that the integral takes the form 



47r" 



1 



ln(a:2 + £2)M2 

a;2(x2+£2) + (x2+£2)2 



Fl(eM)8{x) 



(4.50) 



as £ ^ 0. Any changes in the way the x ** singularities are cut off results only in a change in the function 
F^ieM^. See, for example, (A.9a-A.9b). The first two terms are then regulated in the standard fashion, using 
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(A. 10a - A. 11a). We now insert the results (4.48) and (4.50) in (4.46) and use (4.1) to determine the scale 
dependent part of i^Q (eM) + F^{eM). In this way we obtain the complete limiting form of the cutoff amplitude, 
except for a (5i2(5i3(5i4 term whose unknown numerical constant is called 6o- 

The last graph needed, namely p, is primitively divergent. A special device was used to regulate it in Ref. [1] 
(see below) and it is particularly useful to see that the same renormalized amplitude can be obtained from the 
cutoff procedure. In this discussion below we will use some arguments from the treatment of Ref. [1] which do 
not involve the assumption of formal partial integration. 

The bare amplitude for graph p is 



-L p y-^ii ^ ) 



— 7-J (a;i2,a:;i3,a;i4,£) 

TT 



f {x,y,z,e) = 



1 



1 



1 



1 



1 



1 



g2 y2 



z^ +e^ 



(4.51) 



■^ (x — y) + e^ (y — z) + e^ (z — x) + e^ 
Because the graph is primitively divergent, it is sufficient to cut off only one of the six propagators to obtain 
an amplitude with a well-defined Fourier transform. We therefore add and subtract the product of the last five 
propagators without e and write 



f{x,y,z,s) = 



1 



1 



1 



£2 y 



1 1 

2 z2 (x-y)2 [y^zf ( 



1 



z - 



x)' 



+ r{x,y,z,£) 



(4.52) 



It is easy to show by scaling that the limiting form of the remainder term r(x, y, z, e) is C5{x) 5{y) 6{z) where C 
is given by an integral which is infrared convergent because the integrand is a difference of two terms with the 
same leading infrared behavior. We drop this term henceforth. 

We now write 



f{x,y,z,e) = 



1 



[x"^ + s'^) y'^z'^{x — y) {y — z) (z — x) 



(4.53) 



where we have imitated the special device of Ref. [1] in which the degree of singularity of the first propagator 
was artificially increased. We drop the explicit e^ term in the numerator of (4.53) because it is also shown easily 
to contribute a finite triple S-term as e ^ 0. We now use (2.8) and study 



f{x,y,z,e) 



In [a 



')/e 



y'^z'^ {x — y) (y — z) (z — x) 



(4.54) 



The term in brackets is regular as a; ^ 0, so that the box operator can be integrated by parts without surface 
terms in integrals of (4.54) with smooth test functions such as the Fourier transform studied in Ref. [1]. We 
indicate this partial integration with D, and split the argument of the log to obtain 



,1 ^ 1 • 

f {x,y,z,e} = -^n 



In (a 



e^)M^-lne^M^ 



y'^z^{z-yf {y- zf {z- 



(4.55) 



If we replace In (x^ 



M2 



Xxix^M^ in the first term, the result is just the renormalized amplitude 



of Ref. [1], which was shown there to give ultraviolet convergent integrals with test functions. Note that in 
these integrals the derivatives in D are applied both to the test function and the second factor in (4.55). The 
replacement made above can be justified by studying integrals in which the difference 



In (a 



-e2)M2-lna;2M2 In (l + eVa;^) 



(4.56) 



is inserted in (4.55). Scaling arguments show that such integrals vanish as e ^ 0, because of the structure found 
in Ref. [1] for the y — z subintegrals as a; ^ and because the integrand vanishes faster as a; ^ oo than that of 
(4.55) itself. 
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In the second term of (4.55), we reverse the partial integration and pick up a surface term which is exactly 
that obtained in the rigorous derivation of (2.10) from Green's identity. The second term is then 



-7:'6ix) lim — — -^ -^ ^ = -6^^C(3) 6{x) S{y) 6{z) 

x^u y^z^ (a; — y) (y — z) (z — x) 



(4.57) 



where C,{s) is the Riemann (^-function as found from a study of the limit as a; ^ in Ref. [1]. 

The result of this analysis is given in Table 2 in which we have added —9QTT^g%p5{x)5{y)5{z) to account 
for the finite triple (5-terms dropped above. 



Table 2. The limiting form of hare amplitudes for graphs contributing to the IPI J^-point function. Subscripts 
denote the graphs shown in Fig. 2. The first term in each entry is the renormalized amplitude obtained in 
[1], including the number of permutations of external points Xi required to obtain the full contribution of 
a given graph. Cutoff dependent terms in the limit of small e are then given. The constant B is given in 
(A. 13). Numerical constants bg, bm, ■.-, in 3-loop graphs can be expressed as similar integrals, hut their 
specific form is not required. 



Te{Xi) = -167r^g(5i2(5i3(5i4 



lna;?,M2 



r^(x„£)^ -2.92^12,5340— ^3 



2— perms 



^13 



-247rVlne M (5i2(5i3(5 
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• In^ x"^ Af 2 
r|;(xj,£) -> g^5i25^(x:ii) U f + 2-pemis 



"13 



+ghie^M^T){x,,e) 
+ 127rVln2£2^2^i2<5i3<5i 



r/j(xi, e) 



2g3 



r'^i2<5i3 □ 2 — 



'-34 



d ( I d Inxl^M '^ 
7^ ^ 1 2 2 7^ M 2 

OX^ \ ^lA^li "^^3 ^34 



5— perms 



-2g\TLe^M^V^x^,e 



+ASn^g-' ( - In e^M^ + Ine^Af-" -1-B] (5i2<5i3(5i4 



„4 /, l-pS .^2 71^2 

y\{^^.e) ^ -|-<5i2<534 (n -^ 167r2(l - C(3))<5i3 ) + 2-perms.) 



^13 



35 

2 



3^ 

4 



ln£2M2r^(x„ £) - ^ In^ £^M^T){x,,e) 
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1 d Inxl^M^ 



^13-^14 ^-^3 



^34 
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-Si^n 



"-34 



5— perms 



-TT^g* (Sln^e^M^ + 24 In^ £2^/2 _ ^SBlne^M'^ - 12 - 66^) (5i2(5i3(5: 
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r^(x„e)^f^ <5i2n., / d^^i 



47r2 
1 



4 ln(M-xi)2M' 



2 A/f2 



{u- XlY 



1 a lnx|4M2 



9x3^ I (w — 2^3)2 {u — X4)2 (Ja;^ 2:34 



-(5'^(u-X3) Da 



ln^(7i - X4)2m2 + 2 ln(M - XifM 



2A/f2 



(u — XiY 



-^\ne^M'^Tl{x^,e) + 2g\ne^M^T\{x^,e) 



5— perms 



-.92 ( ^ln2e2M2 + lne2M2 -l-B] T){x,,e) 



2„4 



-247rV 



ln3e2M2 



In e^M^ - Ine^A'r -Bhie-'M^ (5i2(5i3(5i4 



r^(xi,e) ^ ^(^12(534 D 
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ln^a;?3M2 + 21nx?3Af2 



+ 2— perms 
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37r2£^ 
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(Xi - U)2(X3 -m)2 



2— perms 



27r2/ (In^ e^Af^ + 2\ne^M'^ - 2 + k) 512S13S 
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r'' (x- e) ^ -^ 



g_ 

47r4 
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1 (9 Inxl^M^ 



OXn OXo \ XijA OX<2 



^12 



1 d lna;i4M2 

■'-13 '^•^3 •''34 



^ 2. lna;24Af2 9 



1 d lna;§4M2 



^•14 



0X3 \ ^13 0X3 a:;34 



, 2e lnx2 m2 did lna;2 M2 



^34 



0X2 I 2:24 6'a:;2 a;i2 
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47r^ 



<524<5l3 □ 



Iri^xhM'^ + 3\n^xhM^ + 6lnxhM^ 



5— perms 



+glne^M''Tl{x„e) - g^ In e^M^T){x„e) 



J- nl-^i? ^j 



dx'i 



1 9 a; 



^14^34 



X 



14 



d^y 



(y - xiY{y - x^yiy - X4y 



1 d l-a^xl^^M^ + 2liixl4^M^ 



^13-^14 "-^3 



"-34 



,3^2 



2^2 



-y <5i3 D 



In^ xi^M^ + 6 In xf^M^ + 12 lna;f4M^ 



^•14 



11— perms 



+2g\ne^M^rl{xi,s) - 2g^ {ins^M^ + 2lne^M^ -2-2B) T){x^,e) 
16n'^g'^ (In^ e^M^ + 6 In^ e'^M^ + 6(1 - B) Ine^M^ - 12 - 6B + 36„) (5i2<5i3(5i4 



i,, ff'^x X • In^xfgM^ + Sln^xfgM^-einxfgM^ 



n(x„e)-^-^,5i2(534n 



^13 



2— perms 



-gln£2M2r^(x,,e) - y (in^e^^f^ + 21ne2M2 -2-2B) T){x^,e) 
-47r2/ (W e^M^ + SW e^M"^ - Glne^M^ + b^) (5i2(5i3(5i4 



r'(x.,e) 



4g4 ^. Inxf^M^ 

I— I n 



H2 



2 9 2 2 2 

ryt^ ryt^ fvt^ ry*^ rvi^ 

X]^3X]^4X23J-24-i34 



-967r2.g4 (C(3) Ine^M^ + 6^) (5i2<5i3(5i4 



5. Counterterms and Unitarity 

In the previous section, we computed the bare amplitudes for 2- and 4-point functions. These amplitudes 
consist of the renormalized amplitudes of [1] plus terms which diverge as the ultraviolet cutoff e goes to 0. In 
much the same way as in any other regularization method, we now show that the divergent terms can be cancelled 
by local counterterms in the Lagrangian. Specifically, the purely local divergent terms, those involving only delta- 
functions, directly determine counterterms which are added to the Lagrangian. The counterterm vertices in turn 
generate non-local graphs which must cancel the non-local divergent terms in the bare amplitudes. This last 
cancellation is the nontrivial check of the consistency of the entire procedure and the keystone of the proof of 
unitarity. After it is finished we can all go to the beach. 

Of course renormalization conditions are necessarily involved in the process of determining counterterms, 
and we implicitly use the conditions of [1]. The essential role of these conditions is to fix the finite parts of the 
divergent terms, and the consistent cancellation of these by the finite parts of the counterterms below is crucial to 
unitarity. We remind the reader that in the standard form of the renormalization procedure used in this section, 
the coupling g which appears in the Tables is the "physical" coupling within the renormalization scheme of [1]. 
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The simplest illustration of this program is the 1-loop renormalization of the coupling already discussed in 
Section 2. There we saw how the local divergent term in the bare amplitude (2.14) of graph f is cancelled by a 
counterterm in the Lagrangian obtained from the coupling shift (2.15). 

Let us now study the case of the 2-point function. Its total bare amplitude is found adding the contributions 
of graphs a, c, and d from Table 1. The local divergences are cancelled by the following counterterms : 

1. a wave function renormalization, 



-06{x) (Z^ - 1) = -n5{x) 



-^-{l-lne'M-) 



^" '\n'e'M^+lneH4^ + ^ 



2. and a perturbative mass counterterm, 



(5.1) 



Sm Z^5{x) = 5{x) 



1 



1 



;3 



1 



g^TrHns^M^-2Di) 



(5.2) 



Finally, the non-local divergence of graph d, 

Sglne^M'^ tI{x,s) 



(5.3) 



cancels when we take into account in the computation of c, the one loop modification of the coupling (2.15). We 
have, thus, shown how counterterms can eliminate all divergences, local as well as non-local, in the two point 
function. We recover the renormalized amplitude of [1]. 

The same can be done for the four point function at higher loops. The local triple-delta terms are easily 
absorbed into a redefinition of the vertex. One includes, then, all the new counterterm vertices, that is coupling, 
kinetic and mass corrections, in the computation of the amplitude and checks that the non-local divergences 
cancel. Let us stress the non-triviality of such cancellation, which is particularly dramatic for the F^ terms at the 
three loop level. Here, the sum of the 8 3-loop graphs {i — p) contains the term 



4 6 6 ' -' 



(5.4) 



and this is neatly cancelled by the counterterm vertices from both 1- and 2-loop coupling renormalization inserted 
in graph f . In the end, we recover the result of [1]. 

These results can be summarized by adding the counterterms to the Lagrangian (2.1). The result is the 
so-called bare Lagrangian, 



with 



Cbare = -^Z^ {d^^ipf + -Z^Sm^cj)^ + -^ g Z g Z^cj)'^ , 



2 ]\,f2 



ZgZi = 1- -g Ine^lVf + gU - In^ e^M^ + 3 Ine^M^ - 3 - 3B 



16' 
(216B - 44 - 96C(3)) Ins^M^ + c) + O (/) , 



(5.5) 



(5.6) 



Zrf, = 1 - — (1 - Ine^M^) - ^ ( In^ e^M^ + Ine^M^ 



^ )+«(/) 



(5.7) 
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Sni^Z^ = — 



g^- g^{\neH'I^-2D^) 



+ 0{g'). 



(5.8) 



- 7I/2 



We can rewrite 5.5 defining the bare field 0o — Zj (f), the bare coupling qq — gZg and the bare mass jtiq — 5m 
Then, the bare Lagrangian becomes 



1 



1 



-Cbaro = :t (d^M + o'TtO0O + 






-go(t)'< 



• 



(5.9) 



The constant c in (5.6) is the sum of the various numerical constants in the .9^(5i2'5i3(5i4 terms of the 3-loop 
graphs of Table 2. Its value would become relevant only if the calculations of this paper are extended to 4-loop 
order. 

6. Renormalization Group Equations 

The renormalizability of massive (f)"^ theory means that, in any correct regularization procedure with short- 
distance cutoff e and any renormalization scheme with scale parameter M, renormalized and bare IP! n-point 
functions are related by 



ri:l (x„ 5, m, M) = Z^ {go, M, e) T^^^^ (x,, go, mo, e) 



(6.1) 



The right-hand side has a finite limit as £ ^ with physical coupling g and mass m held fixed. 
The fact that 



d 



^^TJ^h'lc i^i^9o, mo, e) 



^jy-barc 







3o,mo,e fixed 



leads directly [4] to the Callan-Symanzik equation 



.^a^ + ^^^^l-^^^-''^-'^^)"^^ 



4")(x„g,m,M) = 



with 



/3(5)=M— 5(go,M,£) 



5o,£ fixed 



7(5) = ^M— lnZ^(go,M,£) 



go -,6 fixed 



^^^"~2^dM'^ (go,mo,M,£) 



(6.2) 

(6.3) 

(6.4) 
(6.5) 
(6.6) 



The relations between bare and physical parameters, such as go = gZg, must be inverted to compute the functions 

Pig), 7(5), Sig). 

The cancellation of divergent terms with counterterms in Section 5 establishes that the differential renormal- 
ization procedure is correct through three-loop order for m = 0, so that (6.1) holds up to three-loop order in the 
massless theory. One can think of the bare amplitudes as defined via (6.1) by substituting g (go), computed from 
(5.6)-(5.7), in the renormalized amplitudes of Ref. [1] and multiplying by Z^" . However, our computational 
procedure also provides quite directly the explicit form of the bare amplitudes. Specifically, the entries in Tables 
1 and 2 are the result of a systematically cutoff computation with Lagrangian coupling g, zero Lagrangian mass 
and unit normalization of propagators. Thus if we simply replace g -^ go in Table 1, we can interpret the entries 
there as rj^"J.^ {xi,go,0,e), and one can verify directly from the table that the apparent M dependence of these 
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amplitudes cancels. To obtain T^J^^{xi,go,ma,e), one simply adds the mass insertion (5.2) or (5.8) rewritten 
in terms of the bare coupling as 

ml = <5m2 - 1 (Igo^ - ^-^g^^ + ©(ff^) . (6.7) 

The sole effect of this is to cancel all quadratically divergent terms in the entries of the Tables, leaving amplitudes 
which clearly satisfy (6.2) because M-^Sm^ = 0. 

The previous arguments establish the validity of the standard formulae (6.4)-(6.5) for the renormalization 
group functions through three-loop order for m = 0. It is a straightforward matter to use (5.6)-(5.7), with proper 
attention to the inversion of go{g) ^^ g{go), and obtain 

(3{g) = 3g' ^g' + (31 + 12C(3)) / + 0(/) (6.8) 



3 

3 

12" " 8' 



7(5) = 49'-|9'VO(/). (6.9) 



Our calculations have probed only the massless theory, and it is clear that the rn-£^ term in (6.3) vanishes as 
771 ^ because there are no infrared divergences. Therefore, we do not discuss d{g) here because it requires 
information about m^Q. The results for /3(g), 7(17) found here in the standard framework of the renormalization 
group equations agree with those of the "experimental" approach taken in Ref. [1]. This provides another check 
that the differential renormalization procedure is correct. 

7. An Alternate Cutoff Method 

The cutoff method used in Sections 1-4, which is based on the damped propagator (1.2), is actually the 
second method we have applied to this problem. In our first approach, which was described briefly in the first 
article of Ref. [1], regularization of bare amplitudes was achieved by the exclusion of small balls of radius e 
about short distance singularities. Integrals involving such cutoff bare amplitudes then converge, and the singular 
contributions as e ^ are quite clearly related to the surface terms dropped in the partial integration rule of 
differential renormalization. 

The systematic rules used in this regularization method were the following: 

1) Each propagator connecting vertices x and y of a diagram is replaced by the cutoff propagator 

where 0(z) is the step function. 

2) For each pair of (internal or external) vertices Xi,Xj not connected by a propagator, the bare amplitude is 
multiplied by an additional cutoff factor (|xi — Xj\ —s). 

Calculations using this approach were generally simpler than in the current method because the differential 
identities of the Appendix of [1] could be used directly. Further, after partial integration, one finds S {\xi — Xj\ ~ e) 
terms when derivatives act on the step function cutoff factors, and these effectively reproduce the surface terms 
which are the crucial issue. Complete results through 3-loop order were obtained, and we found that the cutoff 
amplitudes for each graph could be expressed as the renormalized amplitudes of [1], plus singular terms which 
could be consistently compensated by counterterms in the Lagrangian. The renormalization group functions f3{g) 
and 7(g) were calculated from the cutoff dependence of Zg and Z^, as in (6.4) and (6.5), again with results 
identical to those of [1]. 
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Despite the successful result and relative ease of calculation, we now believe that this method does not 
support the conclusion that the renormalized amplitudes of [1] satisfy perturbative unitarity. To discuss this we 
first note that the cutoff propagator (7.1) has Fourier transform 

A(p,e) = ljo(pe) (7.2) 

where Jo(z) is the Bessel function which is analytic in its argument. Since the only singularity of A{p,e) is the 
standard l/p^ pole, the method would give a plausible argument for unitarity provided that calculations could be 
done using only the propagator cutoff of Rule 1) above. 

In principle, the propagator cutoff is sufficient to make all required integrals converge, but it was technically 
too difficult to do many integrals in this way. Instead we adopted the procedure of performing subintegrals using 
Rule 1), but then substituting the limiting form of this result before studying further singularities of a graph which 
were cutoff by factors from Rule 2). These Rule 2) cutoff factors cannot be described as a modification of the 
Lagrangian which is Hermitean below some cutoff energy scale, and this raises more questions about unitarity. 

In view of the above, one may wonder whether the result of the consistent counterterms mechanism found in 
this method was an accident or whether it encapsulates some truth. We think that the latter is correct, because 
our procedures, albeit somewhat sloppy, were used consistently. Subgraphs of a given graph were handled by the 
same steps as in their initial appearance in lower order. 

Very recently, it has been shown that x-space dimensional regularization can be combined with differential 
identities so as to reproduce several of the lower order amplitudes of [1] plus local counterterms [5]. In higher 
order, this method could lead to a useful relation between the amplitudes of the differential renormalization and 
dimensional regularization procedures. 

8. Concluding Remarks 

We believe that the calculations of Sees. 2-4 have fulfilled their intended goals. Namely, a systematic real 
space cutoff method for 0^ theory has been used to show that through 3-loop order, the bare 2- and 4-point 
correlation functions can be expressed as the sum of the renormalized amplitudes of [1] plus a combination of 
singular and finite terms. This combination can be compensated by adding the traditional counterterms to the 
Lagrangian. Indirectly this demonstrates that the major heuristic rule of the differential renormalization procedure, 
namely formal partial integration, is consistent. Since the cutoff method is based on a damped propagator whose 
Fourier transform (1.4) consists of the usual pole plus a cut whose effects vanish as e^, our results also imply that 
differential renormalization obeys perturbative unitarity. Finally we have shown that the same renormalization 
group functions (3{g) and 7(g) are obtained in the cutoff theory and in the method of [1], and this is an additional 
consistency check. 

An important subsidiary purpose of our work was to convince skeptics that overlap divergences are correctly 
treated in differential renormalization. The results above do demonstrate this since all non-local divergences are 
exactly cancelled by the local counterterms added to the lagrangian. However, it may be useful to restate and 
amplify upon the common belief that overlap divergences are not a problem in real space calculations, because 
subdivergent regions remain distinct and can be regulated before the overall divergence of a graph is studied. Let 
us illustrate this in the case of the most conspicuous overlap graph in our work, the 3-loop cateye graph o. We 
note that the treatment of this graph in [1] started with the expression 

4g4 
^o{xi) = ^Si2634fo{xi3) + 2-perms 

1 /• d*ud^v . ln{u-v)^M^ ^^-^^ 

4: J u'^v'^{x — u)'^{x — v^ (m — w)2 
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This expression can also be obtained from (4.45) by setting e = and using (1.1) to regulate the central bubble 
subgraph. 

There are three subdivergent regions in this graph; namely the 4-dimensional region m ~ w, and the eight 
dimensional regions u ^ v ^ and u ^ v ^ x. The first of these is already regulated by the use of (1.1) which 
is an identity for u — u 7^ 0, and is defined at the singularity by the partial integration rule. This means that the 
divergent 1-bubble subgraph is treated exactly in the same way as the renormalized amplitude for graph f . After 
partial integration ofDu in (8.1) one finds [1] two integrals with S{u) and S{x — u) factors together with a cross 
term, as one can see from (4.46) at e = 0. The cross term is complicated but contains no subdivergences. Indeed 
the S{u) and S{x — u) factors give a partial localization of the 8-dimensional singular regions, and it is not difficult 
to see that the integrand in these regions is treated in [1] exactly as the product of the renormalized amplitude of 
the ice cream cone subgraph (specifically, the first term in the entry for F^ in Table 2 with correct designation of 
variables) times the remaining non-singular propagator factors. The internal integrals d^uct^v are then performed 
leading to the explicit form of fo{x) containing (Inx^M^)"/^^ overall singularities which are easily regulated. 
Independent of a systematic cutoff procedure, our confidence that overlap divergences are properly treated in 
differential renormalization is based on the property that the amplitude in subdivergent regions is regulated in 
the same way as the appropriate subgraph. 
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Appendix 

We include in this Appendix some technical results. Most of them have been used in the calculations of 
several of the graphs discussed in Sections 2-4 of the main text. 

1. The convolution integral 



iJx) = / d^u ^ (A.l) 

is required to evaluate graphs d, g, i, and j . The integral can be done using the fact (see (2.12)) that the Fourier 
transform of 1/(m^ + e^)^ is 2'K^Ko{pe). We can write 



(27r) 



lei^) = / jd^e-^P-An^Kiips) 



-n J d^p—^Kiipe) 
-TT^nl - / dpJi{px)Kl(pe) 



arctanh \ \ 1 



W 



n 2 \ 



-TT^D 



^2 ln2 [(a;2 + 2e2 - x^x^TIe^ /2e^] 
4 a;2 



(A.2) 



This is an exact result. As e ^ we obtain 



I,{x)^-^u'^^^4^. (A.3) 



If we introduce the mass scale M , this can be rewritten as 



■,2 ^2 ]i/r2 ^2 . i„ ^2 ]\/r2 



Inx^M^ tt\ 2,^2 • lna;2Af2 
4 a;2 2 x^ 



Ie{x)^ — ^ D ^^5 + — In e^M^ D 5 + t:'^5{x) In e^M^ . (A.4) 



This is the same result we would have obtained if the limiting form of the regulated bubble 

1 1 . lna;2M2 



D 5 7r2 \ne'M^5{x) , (A. 5) 

were inserted in (A.l), and the integral computed using formal partial integration as in [1]. In other words, it is 
justified in this case to take the limit in the integrand. The reason for this appears to lie in (2.8) in which the 
cutoff bubble amplitude is expressed as Q of a function which has a soft singularity as e ^ 0. The D operator 
can be transferred to the external variables in (A.l), leaving a function which is smooth enough that the e ^- 
limit can be taken in the integrand. This is sufficient for the evaluation of graphs g and i. However, in graphs 
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d and j, where other singular factors multiply I^{x), the more accurate form (A. 2) is required to study the limit 
of the bare amplitudes. 

2. Representations of the distributions S{x) and nS{x) appear throughout our work. The simplest example is 
the basic equation for the cutoff propagator, 

A. ^ 1 1 2e2 
nA(a;,e) = -—^n-r^ ^ = ttt^ ^t^t . (A. 6) 

^ ^ 47r2 (x2+e2) 7^2(^2 +£2)3 ^ '' 

Of course we expect the limiting relation 

2e2 



2 ( ^2 i c-2~\3 



7T^(x^ + S' 



5{x) . (A.7) 



To prove this we integrate this candidate (5-function with a smooth test function f{x) which is damped at long 
distances. We write 



/^V(-)^,(J^ = ^JdSfiey) (^ 



71"^ J (2/2 + 1)3 7r2 y ^ (-y2 _^ ^)3 



It follows from Taylor's theorem that f{sy) — f{0) ^ 0(e) as e ^ so the limit of the second infrared convergent 
integral vanishes. The first integral is easy to evaluate, 



d^y 


^J''l iy^ + ir 
^ , r y'dy 

' h {y' + ir^ 2 


(y2 + 1)3 



Other representations of distributions are 

3e2 



(a;2 + £2)4 2£2 

£2 ln(x2 + £2)M2 



d{x) + —Dd{x) (A. 8a) 



,, ^ „,, ■ 7r^lne'M^ + l)5{x) . (A.8b) 

X''(X^ + £^)^ 

When nS(x) is involved, it is necessary to expand the test function in a Taylor series through second order in 
order to extract the limiting form. Similar relations which hold as £ ^ are 

^ ^ -TT^S{x) (A.9a) 



(a;2 + £2)2 a:;2(a;2+£2) 
ln(a;2 + £2)/£2 ^ ln(x2 + £2)/£2 ^ 



(a;2+£2)2 2,2(2;2+£2) 

ln"(x2+£2)/£2 _ ln"(a;2 + £2)/£ 2 ^^ 
(x2 + £2)2 a;2(a;2 + £2) 



TT-'Six) (A.9b) 

n'^nlSix) . (A.9c) 



3. Differential Regularization identities for cutoff singular functions are useful throughout our calculations. 
These include the following identities, 

' - i^M:f!+£!)Z£! (^,.103) 



(x2+£2)2 4 X 
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(a;2 + £2)3 32 



1 ln(a;2+£2)/e2 

— nn— ^ 77^^ — 



3e2 



(x2 + £2)4 



In(x2 + e2)/e2 ^ W {x^ + e^) j e^ ^ 2\xi{x^ + e^) j e 



(A.lOb) 
(A.lOc) 



As £ ^ 0, we also have 

ln(x2 + £2)M2 1 \x?{x^ + £2)M2 + 2 ln(a;2 + £2)M2 7r2 2 2 2 



a;2(x2 + £2) 



(in^ e^M'^ - 2) (5(x) 



In^(x2+£2)M2 _\_ ln^(a;2 + £2)Af 2 + 3 ln^(-^2 _^ ^2)^^-2 ^ g ^^(^2 _^ g2)jy^ 

x2(a;2+£2) ' "12'-' ^^2 



-^2(lln3^2^^2_l')^(^)_ 



(A.Ua) 



(A.llb) 



4. Triangular structures have the schematic form, 

[a representation of b{x — j/)] > 



smootfi /(j/) 



(j/2 + £2) or J/2 



smooth (/(x) 



(a;2 + £2) or x"^ 



Such structures first appear in the 2-loop graph h, and we also encounter them in many 3-loop graphs. The basic 
strategy to obtain their limiting form, which involves integration with a test function f{x,y) of two variables was 
discussed in connection with graph h. Using this strategy one obtains the limiting forms 



1 



1 1 

a;2 + £2 y2 



1 



(x — y)2 + £2 

1 ln(y2+e2)M2 



— 47r^(5(x — y) 



1 



2/2 (y2 ^ g2) 



- 47r4(5(x)(5(2/) 



(A.12a) 



D 



(x — y)2 + £2 
1 

(x — J/)2 + £2 
1 



-47r^(5(x — y) 



ln(2/2 + £2)M2 



(x — t/)2 + £2 



1 ln"(y2 + e2)/£ 



ln(x2+£2)/£2 1 

j,2 _|_ g2 2/2 



-47r (5(x — y) 



2/2 (j/2 _|_ g2) 

ln"(y2+£2)/£2 
2/2(y2_^£2) 



+ 47r*(ln£2M2 - B)(5(x)(5(j/) 



(A.12b) 
47r^BJix)6iy) (A.12c) 



2e/ , lny2 4.e2 /£2 

-47r 5 X - y) — „, „ „, + 47r (5(x)(5 y) 



(A.12d) 



The coefficient of 5{x)S{y) in (A. 12b) is simply the double integral of the last term in (3.16). This leads to 
the following expression for B, 



B^^Jd'xd'yU 



1 



(X - J/)2 + 1 
1 



1 



ln(y2 + 1) 



x2(x2 + 1) 2/2 



= _1 /'^4 .4 1 ^__ 

TT^J \x-yf + l X2(x2 + 1) 



y2 + l 



9 3^ V3 



-1.171953 



(A.13) 



One can go from the first to the second line of (A.13) by integrating D by parts back onto the log. Finally, this 
last integral is computed using the standard Feynman parametrisation. Similar expressions can be found for the 
numerical constants i3„. For example, the coefficients of the triple delta terms can be read off by integrating the 
difference of both sides of the equalities (A. 12). 
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